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Abstract. The Cauchy problem for the Fokker-Plank-Kolmogorov equation with a nonlo- 
. cal nonlinear drift term is reduced to a similar problem for the correspondent linear equa- 

04 ] tion. The relation between symmetry operators of the linear and nonlinear Fokkcr-Plank- 

Kolmogorov equations is considered. Illustrative examples of the one-dimensional symmetry 
operators are presented. 
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g '. 1 Introduction 

By definition, symmetry operators leave invariant the solution set of an equation and allow 
^ . to generate new solutions from the known ones (see, e.g. [H [2]). The finding of symmetry 

operators is an important problem, but it rarely can be solved explicitly because the equations 
that determine symmetry operators are complicated and nonlinear. Therefore, special types of 
symmetry operators are of interest. The most affective approach is developed in the framework 
of the group analysis of differential equations [31 El El [6], [7] where the Lie groups of symmetry 
operators are considered. The Lie group generators (related to symmetries) are obtained from 
the determining linear equations which can be solved in a regular way when the symmetries are 
differential operators. The symmetries of differential equations can be considered in the context 
of differential geometry [U [9] . 

The calculation of symmetries for integro-differential equations is a more complex problem 
because there are no general way of choosing an appropriate structure of symmetries. The 
finding of symmetry operators for a nonlocal equation is usually a hopeless task. Under these 
circumstances, examples of symmetry operators for a nonlinear equation with nonlocal terms 
are of mathematical interest. 

In this work we consider an approach that can be used to obtain such examples for the 
Fokker-Planck-Kolmogorov equation (FPKE) of special form with a quadratic nonlocal nonli- 
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nearity 

{-dt + eA + dis(y(x,t) + xJ W(x,y,t)u(y,t)dy^u(x,t) = 0, (1.1) 

where 

V(x, t) = Ktx, W(x, y, t) = K 2 x + K 3 y. (1.2) 

Here, t € R , x = (xi, . . . , x n ) T G M n , j/ = (yi,...,y n ) T S M n are independent variables; 
(xi, . . . , x n ) T means a transpose to a vector or a matrix; (fx* = dx\ ■ ■ ■ dx n ; the dependent variable 
u(x,t) is a real smooth function decreasing as ||x| — > oo; K±, K2, K% are arbitrary constant 
matrices of order n x n; e and >r are real parameters; dt = d/dt; dg = d/dx is a gradient operator 

n 

with respect to x; A = <9f<9j = d 2 jdx\ is a Laplace operator. 

i= l 

The operator of equation (jl.ip is quadratic in independent variables and in derivatives and 
it has a nonlocal quadratic nonlinear term. 

This equation serves as a simple example of a class of "near- linear" nonlocal equations |lUj . 
such that they are nonlinear but the integrability problem for them can be reduced to seeking 
a solution of appropriate linear equations. Nonlinear equations of such type regularly depend 
on the nonlinearity parameter and they possess solutions which go into solutions of the linear 
equation as the nonlinearity parameter tends to zero. 

Equation (jl.ip arises in mathematical problems and it can be used in physical applications. 
In particular, the FPKE (|l.lj) . (jl.2p describes the leading term of the asymptotic solutions 
constructed in [11] in the framework of the formalism of semiclassical asymptotics |12j [T3] for 
equation (jl.ip . in which V(x,t) and W(x,t) are real vector functions of general form. 

The semiclassical approximation is widely used in nonlinear mathematical physics, providing 
a possibility of constructing explicit asymptotic solutions for mathematical physics equations 
coefficients of which can be arbitrary smooth functions and derivatives of dependent variables 
are assumed small. Most of the equations solved by semiclassical methods are not exactly 
integrable. For these equations, semiclassical methods offer a unique opportunity to investigate 
them analytically. 

A method of semiclassical asymptotics based on the formalism of the Maslov complex germ 
[T2j fT3j HH [151 [E] has been developed for a many-dimensional nonstationary Hartree type equa- 
tion with nonlocal nonlinearity in a class of functions localized in a neighborhood of some phase 
curve [17j [TBI fT9[ [20] . This class of functions has been called the class of trajectory-concentrated 
functions (TCF). The Hartree type equation whose operator is quadratic in independent vari- 
ables and derivatives provides another example of the class of "near-linear" nonlinear equations 
similar to the FPKE ([TTjl . 

The symmetry analysis area may be augmented by the study of the symmetry features of 
the semiclassical approximation as the semiclassical methods are hoped to result in a new kind 
of symmetries for mathematical physics equations. The group properties of the semiclassical 
approximation were considered in quantum mechanics and in some models of the quantum field 
theory [21] . The semiclassical method for solving the Cauchy problem in the class of TCF's has 
been developed for the Hartree type equation [T71 [181 121 ED] and for the one-dimensional FPKE 
[111 [22j 123] . For the Hatree type equation, symmetry operators have been found in the TCF 
class. 

A nonlinear FPKE was used to analyze stochastic processes in various physical phenomena. 
In this connection, the following works where the Fokker-Plank-Kolmogorov equation with 
the nonlinear drift term similar to that in (jl.ip was considered deserve mention (see also [10| 
and references therein). M. Shiino and K. Yoshida studied noise effects and phase transitions 
effects involving chaos- nonchaos bifurcations [21] in the framework of nonlinear Fokker-Plank 
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equations. These equations are shown to exhibit the property of H theorem with a Lyapunov 
functional that takes the form of free energy involving generalized entropies of Tsallis |25j . 
In [26] the stochastic resonance phenomenon is discussed; in [27] binary branching and dying 
processes were studied. The evolution of quantum systems was described by means of nonlinear 
FPKE's [28] where the nonlinearity reflects the quantum constraints imposed by the Bose and 
Fermi statistics. 

The paper is organized as follows. In Section [2] the nonlinear FPKE is presented with ne- 
cessary notations and definitions. The Cauchy problem for the nonlinear FPKE is reduced 
to a similar problem for the corresponding linear FPKE in the class of functions decreasing at 
infinity via the Cauchy problem for the first moment vector of a solution of the nonlinear FPKE. 
With the help of the Cauchy problem solution we construct a nonlinear evolution operator and 
the corresponding left inverse operator in explicit form for the nonlinear FPKE. 

In Section [3] a general class of nonlinear symmetry operators is considered for the nonli- 
near FPKE. The symmetry operators are introduced in different ways, in particular using the 
evolution operator and the left inverse operator. Examples of one-dimensional symmetry opera- 
tors are given in explicit form as an illustration. In Conclusion the results are discussed in the 
framework of symmetry analysis. 



2 The Cauchy problem and the evolution operator 

In our consideration, the key part is played by the Cauchy problem for the FPKE (jl.ip . (jl.2p 
in the class of functions u(x,t) decreasing as \\x\\ — > oo at every point of time t ^ 0. To be 
specific, we assume that u(x, t) belongs to the Schwartz space S in the variable x € M n and 
regularly depends on t, i.e. u(x, t) is expanded as a power series in t about t = 0. Obviously, 
equation (II. lh can be written in the form of the balance equation 

dtu(x, t) = d$B(x, t, u), 

where 

B(x, t, u) = edgu(x, t) + V(x, t)u(x, t) + x / W(x, y, t)u(y, t)dyu(x, t). 

Then according to the divergence theorem, we obtain that the integral J* K?l u(x, t)dx conserves 
in time t for every solution u(x,t) of equation (jl.ip . Therefore, taking the initial function 
u(x,0) = j(x) to be normalized, L n j(x)dx = 1, we can assume 

/ u(x,t)dx = l, t^O, (2.1) 
Jm. n 

without loss of generality. We do not pay special attention to the positive definiteness of solutions 
of the FPKE (jl.ip . leaving this requirement for specific examples (see |29j for details). 
Let us write equations (jl.ip . (jl.2p in equivalent form 

{-d t + H nl {x, t; X u (t))}u(x, t) = 0, (2.2) 
where the operator H n \ reads 

H nl (£,t;X u (t)) =eA + d g (Ax + xK 3 X u (t)), (2.3) 
the matrix A is 



A = K x + xK 2 , 
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and the vector 

X u (t) = I xu(x,t)dx (2.4) 



is the first moment of the function u(x,t). With the obvious notation X u (t) = dX u (t)/dt, we 
obtain immediately from (pl2" ]H |ll4" ]) . and ((2TT]) 

X u (t) = -(A + xK 3 )X u (t). (2.5) 

Equation (|2.5p can be considered the first equation of the Einstein-Ehrenfest system (EES) that 
describes the evolution of the moments and centered high-order moments of a solution u(x, t) of 
the FPKE (jl.ip with the vector functions V(x, t) and W(x, t) of general form. The total EES 
for moments of all orders was derived in constructing approximate semiclassical solutions for a 
one-dimensional FPKE in [IT] . 

2.1 Solution of the Cauchy problem 

Let us set the Cauchy problem for equation (|2.2|) : 



u(x,0) = 7(x), j(x) € S, / ^{x)dx = 1. (2-6) 

JR™ 

Then we have the induced Cauchy problem for the vector X u (t) 

X u (0) = X J= xj(x)dx (2.7) 
determined by (|2.5[) . 

The nonlinear Cauchy problem (|2.2p . (|2.6p is reduced to a linear one as follows. For a given 
initial function j(x) (|2.6|) . we can seek a solution of the Cauchy problem (|2.5|) . (|2.7|) indepen- 
dently of the solution of equation (|2.2p and obtain the vector X u (t) having the form of (|2.4p 
due to the uniqueness of the Cauchy problem solution. Let us introduce a function w(x, t) by 
the equality 

u(x,t) = w(x - X u (t),t). (2.8) 

By substitution of (|2.8p in (|2.2p we obtain for the function w(x,t) a linear equation: 

-d t w(x,t) + Lw(x,t) = 0, (2.9) 
L = eA + drgAx. (2.10) 

From and $L8\i we have 

w (x,0) = 7(f) = 7(x + X 7 ) (2.11) 

and 

/ 7(f)df=l. (2.12) 

JR™ 

Equation (|2.7|) results in 

/ xj(x)dx = 0, (2.13) 
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i.e. the function j(x) is centered. Obviously, the integral J Rn w(x,t)dx conserves in time t for 
any solution w(x,t) of equation (|2.9|) ; then from (|2.12|) we have 



w(y,t)dy = 1. 



Therefore, the nonlinear Cauchy problem (|2.2p . (|2.3|) . (|2.6p can be solved as follows. First, 
for a given initial function j(x) (|2.6h we solve the linear Cauchy problem f)2.9h . (12.10p . (12.111) 
with the initial function 7(x) normalized by condition (|2.12p and centered by (|2.13p . Second, 
we find the vector X u (t) by solving the Cauchy problem (12. 5p . (12. 7p . Then, the solution of the 
nonlinear Cauchy problem ([22D, pT3]) . (j2"U|) is given by ([2T8]) . 

Equation (|2.9p is known (see, e.g., [2]) to have a solution in the form of a Gaussian wave 
packet: 



w(x, t) 



'detQ(t) 
(2vre) n 



exp 



— £ T Q(i)£ 



(2.14) 



where Q(t) is a symmetric positive-definite matrix of order n x n. Substituting (|2.14|) in (|2.9|) . 
we obtain 



d 



x T Q(t)x + 2 x J (Q(t)) z x - x T A T Q(t)x - x T Q(t)Ax - e— log det Q(t) 

+ 2eTr (- Q(t) + A) = 0. 

Here Tr A is the trace of the matrix A. Equating the coefficients of equal powers of x, we have 

Q(t) + 2(Q{t)) 2 - A J Q(t) - Q(t)A = 0, (2.15) 
d 



dt 



log det Q(t) + 2 Tr (- Q(t) + A) = 0. 



Let us take Q(t) in the form 
Q(t)=B(t)(C(t)y 1 1 



(2.16) 



where B(t) and C(t) are matrices of order n x n. On substitution of (|2.16p in (|2.15p we can 
write 



B{t) = A J B(t), B(0) = B , 

C(t) = 2B(t) - AC(t), C(0) = C , 



(2.17) 



where Bq and Co are arbitrary constant matrices of n x n order. We call equations (|2.17|) 
a system in variations in matrix form. 

For the one-dimensional case, the linear equation (12. 9p takes the form 



{-d t + edl + d x Ax}w(x, t) = 0, 

where d x = d/dx. 

The solution (|2.14p of equation (|2,18p reads 



(2.18) 



w(x, t) 



1 B(t) 
2ireC(t) 



cxp 



Bit) 
2eC{t)' 



where B(t) and C(t) are a solution of the system in variations (|2.17p in the one-dimensional 
case. For t = we have 



w(x, 0) = 7(0;) 



Bn 



2vreC 



exp 



B t 



2eQ 



2 

X 



(2.19) 
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Notice that the function 7(2;) (|2.19p is normalized and centered: 

-+00 



; y(x)dx = 1, / x^{x)dx = 0. 
Then the function 



u(x, t) = w(x — X u (t),t) = a ^ exp 



B{t) {x-X u {t)f 



(2.20) 



2ireC(t) ^ [ 2eC{t) 
will be a solution of the nonlinear equation 

{-d t + ed 2 x + d x Ax + xK 3 X u (t)d x }u(x, t) = 0, (2.21) 
with the initial condition 

B 



[/ 2^eC Q 6XP 



2eC 



(x - x^y 



(2.22) 



The vector X u (t) = £ 

-00 xu ( x i t)dx in equation (12.2ip satisfies the condition 
X u (t) = -(A + *K 3 )X u (t), X u (0) = X y . (2.23) 

2.2 The evolution operator for a nonlinear FPKE 

Let us rewrite the solution of the above nonlinear Cauchy problem in terms of the corresponding 
nonlinear evolution operator. 

Let G\i n (t, s, x,y) be the Green function of the linear equation (|2.9p . i.e. 



w(x,t) = G\ in (t,s,x,y)j(y)dy. 

Substituting x for x — X u (t), we find 

w(x- X u (t),t) = I Gi in (t,s,x-X u (t),y)j(y + X 7 )dy. 
According to (|2.8p and (|2.1ip . the function 

u(x,t)= G n i(t,s,x,y,j)j(y)dy= G\ in (t, s,x - X u (t),y - X^)>y(y)dy (2.24) 



is a solution of the nonlinear equation (|2.2p . (|2.3p with the initial condition (|2.6p . Therefore, 

G n i(t,s,x,y,~/) = Gn n (t, s,x - X u (t),y - X 7 ) (2.25) 

is the kernel of the evolution operator for the nonlinear equation (|2.2p . Here j(x) = u(x, s) and 
the initial time t = is replaced by t = s. 

Suppose that a solution of the system in variations (|2.17p has the form 

B{t) = Mi(t, s)B , C{t) = M 2 (t, s)B + M 3 (i, s)C , (2.26) 

where M\(t,s), M2(t,s), and M 3 (t,s) are the matrix blocks of order nxnof the matriciant 
(evolution matrix) of the system in variations (|2.17|) . These matrix blocks must satisfy the 
condition 

M = AM, M(a) = I 2 nx2n, (2.27) 
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where 

-Kit s fM x (t,s) \ / At \ 

M = M C' S >=U(M) M 3 (M)J' A= U»<„ -aJ' 

Set the Cauchy problem for |[2.26]) : 

B(s) = Bl = B , C(s) = C = 0. 
The Green's function of equation (12, 9p is known and can be taken as (see, e.g., [2]): 

Gvm(t, s, x, y) = . 

^/(27r e )"det[M 2 (i, S )(M 1 (i, S )) *] 

x exp j-^-(x - M 3 (t, s)y) T [M!(t, a)(M 2 (*, s))" 1 ]^ - M 3 (t, . 
Then the nonlinear evolution U(t, s, •) operator (|2.24p reads 

U(t,s, 1 )(x)= f G nl (t, S ,f,y, 7 ) 7 (y)(iy. (2.28) 
The left-inverse operator C/~ 1 (t,s,-) for the operator (|2,28p is 

U~ 1 {t,s,u){x)= I G~^(t,s,x,y,u)u(x,t)<M. (2.29) 



Here (t, s , x, y, u) is the kernel of the left-inverse operator, which is obtained from (|2.25p if 
we substitute t for s and s for i. The explicit form of this function is 



(2^e)"det [M 2 (s,t)(M 1 {s,t)y 1 ] 
x exp{ - l(x-X 7 -M 3 (s,t)(y-X u (t))) T [M 1 ( S ,t)(M 2 ( S ,t))- 1 ] 
x (x-X 7 -M s (s,t)(y-X u (t)))Y 
The Green's function for the one-dimensional equation (|2.18p is 

1 Mi(i,s), 



Ci in (M,x,y) = W 27rgM2(M) exp 



2eM 2 (t,s) 



(aj-M 3 (t,s)y)' 



where Mi(i, s), M 2 (i, s), M 3 (t, s) are solution of the system in variation (12.27P in the one- 
dimensional case. 

The kernel of the evolution operator for the nonlinear equation (|2.2ip takes the form 

Gni(t,s,x,y,-f) = G hn {t,s,x - X u (t),y - X y ) (2.30) 



1 Aft(t,s) 
27reM 2 (t,s) 6XP 



S \( x -X u (t)-M 3 (t,s)(y-Xj) 2 



2eM 2 (t,s) 



where X u (t) satisfies equation (12.231) . The evolution operator (|2.28p with the kernel (|2.30p is 
written as 



u(x, t) = U(t, s, 7)(z) = / GW(i, s,x- X u (t),y - X^(y)dy. (2.31) 

J — oo 

Here the function u(x,t) having the form of (|2.3ip is a solution of equation (|2,21l) . 
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Notice that direct calculation of the action of the evolution operator (|2.28p on the func- 
tion 7(2/) having the form of (|2.22|) gives the function (|2.2U|) : 



U(t,s, 7 )(x) 



+00 



x exp 



2ireM 2 (t,s) 
1 Mi(t,s) 



2eM 2 (t,s) 



( x - X u {t) - Mz{t,s){y - X^))' 



2neC 

I B(t) 
2neC{t) 



exp 



Bn 



exp 



2eC 
Bit) 



(y - x,f 



dy 



2eC{t) 



(x - X u (t)f 



Conversely, the action of the left-inverse operator (|2.29p on the function (|2.20p in the one- 
dimensional case gives the function (|2.22j) : 



U~ l (t,s,u)(x) 



+00 



x exp 



1 MiM) 
2ireM 2 (s,t) 

1 Mi(s,t) 
2eM 2 {s,t) 



{x - X 7 - M 3 (s,t)(y - X u (t))y 



1 B(t) 
2ireC(t) 



exp 



B(t) 
2eC{t) 



(y - X u (t)f 



dy 



2neC ( 



exp 



go 
'2eC 



(X — X-y 



Because the solution of the nonlinear equation (|2.2p is reduced to seeking the solution of linear 
equation (|2.9p in terms of the moment X u (t) (12. 5p . the symmetry operators of these two equations 
are closely connected. 

Equation (|2.9p with the operator L having the form of (I2.10|) is a special case of the linear 
evolution equation quadratic in derivatives d$ and independent variables x. This equation is 
known to be integrated in explicit form (see, e.g., [2]) which in turn leads to integrability of the 
nonlinear FPKE (|2.2I) . (12. 3p according to (j2.8f) or (|2.24p . The basis of solutions and the Green's 
function for equation (|2.9p can be constructed with the help of symmetry operators of special 
form following, for example, [2j [30l IHT]. 

Consider the symmetry operators for (|2.9p and (|2.2p . 



3 The symmetry operators 

The symmetry operators for equation (jl.ip can be found in various ways following the general 
ideas of symmetry analysis [21 HI O [6] . 

3.1 The determining equation and intertwining 

Let us construct for a function 7(3:) of the space S the function u(x, t) of (j2.8f) using the solutions 
of the Cauchy problems (HTTP for the vector X u (t) and pity . OTTTD for w(x,t). 

Let us take an operator a(x) acting in the space S as the initial operator for a time depending 
operator A(x, t): 



A(x,0) = a(x). 
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The function 

ja(x) = — a(x)j(x), Ola = [ a(x)j(x)dx. (3-1) 
"A JR" 

determines the vector X^ A by formula (|2.7p . where j(x) is replaced by 7^4 (x). Taking the 
vector X lA as the initial condition for equation (|2.5p . we find the vector X UA {t). 
Obviously, the operator A(x,t) determined by the conditions 

( - d t + Hni(x, t, X UA (t)))A(x, t) = B(x,t)(-d t + H nl (x, t, X u (t))) , (3.2) 
A(x,0) = a(x) (3.3) 

is a symmetry operator for equation f|2.2|) . Here B(x,t) is an operator such that B(x,t)(0) = 0. 
This operator plays the part of a Lagrangian multiplier and it is determined together with A(x, t). 

Equation (|3.2|) is the determining equation for the symmetry operators of equation (|2.2p . In 
general, (|3.2p is a nonlinear operator equation. But, given the initial function 7(2?) and the initial 
operator a(x), we can find the vectors X u (t) and X UA {t) without solving the equation (|2.2p . On 
substitution of these vectors in f)3.2|) the operators H n \(x,t,X u (t)) and H n \(x,t,X UA (t)) become 
linear. Then we can assume that a(x), A(x,t), and B(x,t) are linear operators in (|3.2p . 

Notice that in the general case the operators A and B depend on X u (t), X UA (t), i.e. 

A = A{x,t;X u (t),X UA {t)), B = B(x, t; X u (t), X UA (t)). 

If B(x,t) = A(x,t) then A(x,t) is called the intertwining operator for the linear operators, 
satisfying the condition 

{-dt + H nl (x, t, X UA (t)))A(x, t) = A{x, t)(-dt + H nl {x, t, X u {t))) , (3.4) 
A(x,0) = a(x). (3.5) 

If the operator a(x) is given in (|3.2p . (|3.3p (or in (|3.4p . (|3.5p ). these conditions are the Cauchy 
problems determining the operator A(x,t). 

Now, let u(x,t) be a solution of the Cauchy problem (|2.2p . (|2.3p . ()2.6j) . and the operator A 
is determined by the solution of the Cauchy problem (|3.2p . (|3.3p or (|3.4|) . (13. 5p . Then we 
immediately obtain that the function 

ii A (i,t) = 4(x,t;i u (i),4(t))^,t) (3.6) 
is a solution of the Cauchy problem 

d t UA(x,t) = H n i(x,t,X UA (t))uA(x,t), 
u A (x,0) = ta(x) = — S(x)7(x). 

Therefore, the operator A(x, t; X u (t), X UA (t)) is a symmetry operator of the nonlinear FPKE 

(22]), (J23D. 

Notice that the operator A(x, t; X u (t) , X UA (t)) in (]3.6p is nonlinear due to the presence of the 
vectors X u (t), and X„ A (i). 
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3.2 Symmetry operators of nonlinear and linear FPKE 

We now deduce a relation connecting two solutions of the nonlinear FPKE (|2.2p . (|2.3p using 
a symmetry operator of the linear equation (|2.9f) . (|2.10p . This relation can be considered a sym- 
metry operator for the nonlinear FPKE. 

To this end consider equation (|2.8j) which connects the nonlinear Cauchy problem (12, 6p with 
the linear Cauchy problem (|2.1ip . 

Let A(x,t) be a symmetry operator of the linear equation (|2.9p . Then the function 

wa(x, t) = -^—A(x, t)w(x, t), ola = f A(x,0)w(x,0)dx, 
ola Jm™ 

is another solution of the linear equation, which is determined by (|2.9f) . (12.101) . For t = we 

have 

w A (x,0) = — A(x,0)7(») = ja(x) 

OLA 

and 



/ jA(x)dx = 1 
Jw 

which leads to normalization of the function wa (x, t) : 
wa(x, t)dx = 1. 



On the other hand, the function wa(x, t) is not centered for a symmetry operator A(x, t) of 
general form. In other words, for t = the vector 

Xa = xwa(x, 0)dx = / xjA(x)dx 

is nonzero. 

To construct a solution of the nonlinear equation (|2.2p , which would correspond to the solu- 
tion WA(x,t) of the linear equation (|2.9p with the use of relation (|2,8p . the function WA(x,t), 
being a solution of equation (|2.9p . should be centered. 

We can immediately check that equation (|2.9p is invariant under the change of variables 

t' = t, x' = x — l(t), where l(t) satisfies the condition l(t) = —Al(t). 
Taking into account this property, let us introduce a vector 



1-A(t) = / xwA(x,t)dx 
which satisfies the Cauchy problem 

t A {t) = -Al A (t), f A (0) = X A . 
Then the function 

w A {x, t) = w A {x + I a (t) , t) 
satisfies equation (|2.9p and the initial condition 

wa{x,0) = w A (x + Aa,0) = j A {x + X A ). 
The function 7a(£ + Xa) is normalized and centered. The same is true for WA(x,t). 
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Following (|2.8p , we now construct a solution va (x, t) of the nonlinear FPKE (|2,2p related to 
w A (x,t). Consider a vector Y(t) such that 

Y(t) = -(A + xK 3 )Y(t), F(0) = A A . 

Immediate check shows that the function 

v A (x, t) = w A (x - Y(t),t) 

satisfies the equation 

{ - d t + eA + d^(Ax + xK 3 Y(t))}v A (x,t) = 0, 
v A (x, 0) = w A (x - Xa, 0) = ia{x). 

Notice that 

X VA (t) = / v A (x,t)xdx, X VA (0) = / xjA(x)dx = X A , 
X VA (t) = -(A + xK 3 )X VA (t), 

then 

Y(t) = X VA (t). 

Therefore, va(x, t) satisfies the nonlinear FPKE (|2.2p . The relation between the solutions u{x, t) 
and va(x, t) reads 

v A (x, t) = w A {x - Y(t),t) = w A (x - Y(t) + l A (t),t) 





= ^-A(x- 

OA 


Y{t) + l A {t),t)w{x - Y(t) + T A (t),t) 






Y(t) + f A (t),t)u{x - Y(t) + T A {t) + ^„(t), t). 


This equation 


determines 


a symmetry operator A n \ of the nonlinear FPKE (|2.2 


ua(x, t) 


= v A (x,t) = 


E 2 n i(f,t)u(f,t) 




= —A(x- 

OLA 


+ f A (t),t) u {x - Y(t) + T A {t) + ^„(t), t). 



(3.7) 

3.3 Symmetry operators in terms of an operator Cauchy problem 

Let us reformulate the construction of symmetry operators for the nonlinear FPKE (|2.2p . (|2.3p 
in terms of an operator Cauchy problem. 

Consider the nonlinear Cauchy problem (|2.2p . (|2,6p associated with the Cauchy problem 
([23j) . ([277] for the vector X u (£) having the form of (|2T3]l . 

With an operator 

a(x) :S^S (3.8) 

acting on the initial function 7(2) G 5 of the Cauchy problem (12. 6p . we define a function 7^4 (x) 
of the form (|3.ip . which is taken as an initial condition for the Cauchy problem for a func- 
tion u A (x, t) 

{ - d t + H nl {x,t,X UA {t))}u A {x,t) = 0, 
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u A (x,0) = j A (x), 

where the vector X UA (t) is determined by the conditions 

4(t) = -(A + xK 3 )4(i), 

X UA (Q) = X 1A , X 1A = I xj A (x)dx. 



Notice that given the function 7(3;) and the operator a(x), we can find X UA (t) not finding 
a solution of the FPKE (p^j) . 

We can immediately verify that the two functions 

w{x,t) =u{x + X u (t),t), (3.9) 
w A (x,t) = u A {x + X UA (t),t) (3.10) 



are solutions of the linear equation (|2.9p and the initial conditions are 
w(x, 0) = 7(x + Xy), w A (x, 0) = 7a(x + X 7A ). 

Define a linear operator A(if, t) by an operator equation 

[-d t + L(x,t),^(x,t)] =0 
with the initial condition 

A(x,0) = a(x), 
where L(x,t) is defined in (|2.10p . It can be shown that 

A(x,t)w(x,t) = A(x + l A (t),t)w(x + l A (i),t). 
Due to the uniqueness of the Cauchy problem solution, we have 

w A (x,t) = A(x,t)w(x,t). 
In view of <^M), (l3J0|) we have 

u A {x + X UA (t),t) = f (x, t)u(x + X u (t),t) 



or 



u A (x,t) = A(x - X UA {t),t)u(x - X UA (t) + X u {t),t). (3.11) 
This relation defines a symmetry operator A n \(x, t) of the nonlinear FPKE (|2.2p : 
u A (x , t) = A^x, t)u(x, t). 
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3.4 Symmetry operators in terms of an evolution operator 

Using the evolution operator (|2.28p and left-inverse operator (|2.29[) . we can obtain symmetry 
operators for the nonlinear FPKE (|2.2p . 

Let a(x) be an operator (|3.8[) acting on an initial function j(x), and u(x,t) is the solution of 
the Cauchy problem (|2.2p . (12, 6ft . Then the function 

u A (x,t) = U(t, s,aZ7 _1 (i, s,u))(x) (3.12) 



is a solution of the nonlinear FPKE corresponding the initial function ja(x) of the form (|3.ip . 
Equation ()3. 12j) defines a symmetry operator A n \ for the nonlinear FPKE (|2.2p : 

ua(x, t) = A n \(x, t)u(x, t). (3.13) 
The one-dimensional case of ()3. 12[) reads 

UA{x,t) = A n \(x,t)u(x,t) = U(t, s,aC/~ 1 (t, s,u))(x), (3-14) 

where U(t, s, •) and U' 1 ^, s, •) are determined by (12T25D and (12T2H in the one-dimensional case. 
Consider an operator a(x, t) of the form 

a(x,t) = M 1 (t,s){x-X u {t)) + (eM 2 {t,s)+M 3 (t,s))d x , 
a(x, s) = x — X-y + d x , X~f = X u (s). 

Here Mi(t,s), M2(t,s), and M 3 (t,s) are solutions of the system in variations (|2.27p in the 
one-dimensional case. Then for (|3.14p we have 

ua(x, t) = lim U(t, s, [dg + z — X~ f ]U~ 1 (r, s, u))(x) 

T — >t 

/+oo r+oo 
dy / d2G n i(t,s,x,z,7A)[9 z + z-X 7 ]G~ 1 (T,s,z,y,u)u(y,T) 



OO 



r^t 2vre y M 2 (t, s)M 2 (s, r) 
x exp { - l(x - X UA (t) - M 3 (i, s)(z - X iA)) 2 ^^}[d z + z - X y ] 

x exp { - 1 (z - X, - M 3 (s, r)(y - X u (r))) 2 ^^}u(y, r) 

= [Mi(i, s)(x - X UA (t) + M 3 (t, s)X lA ) - X 7 + (eM 2 {t, s) + Af 3 (t, s))d x ] 

x u(x + X u (t) - X UA (t) + Af 3 (t, s)(X 7A - X 7 ),t) 
= a(x + X u (t) - X UA (t) + M 3 (i, s)(X jA - X^),t)u(x + X u (t) - X UA (t) 

+ M 3 (t,s){X JA -X 7 ),t). 

Therefore, we have 

u A (x,t) = A nl (x,t)u{x,t) =a(x + X u (t)-X UA (t) + M 3 (t,s)(X yA - X^),t) 

xu(x + X u {t)-X UA (t) + M 3 {t,s){X JA -X y ),t). (3.15) 

In calculating the symmetry operators we have used the following relations: 

M 1 (t,s)M 1 {s,T) =M 1 (t,r), M 2 (t,r) = M 1 (s,T)M 2 (t,s) + M 2 (s,r)M 3 (t,s), 

which follow from (|2.27[) in the one-dimensional case. 
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Equality (|3. 15 j) determines a symmetry operator in explicit form for equation (|2.2ip . This 
symmetry operator generates solutions of equation (|2.21|) . Let us illustrate this by an example. 
By acting with the operator (|3.15j) on the function (12.20|) . we obtain 

u A (x,t) =a(x + X u (t) - X UA (t) + M 3 (t, s)(X 1A - X^),t)x 
x u{x + X u (t) - X UA (t) + M 3 (t, s){X JA - X 7 ),t) 

= W)( CQ ~ (X " + M3(t ' S)iXlA ~ Xl)) 

x u{x + X u (t) - X UA (t) + M 3 (t, s)(X jA - X 7 ), t). (3.16) 

The function UA(x,t) (|3.16p is a solution of the equation 

{-d t + ed 2 x + d x Ax + xK 3 X UA (t)d x }u A {x, t) = 0. 

Notice that the symmetry operators determined by (|3.13p are consistent with the operators A3. Tf) 
and (|3.1ip in the one-dimensional case. Moreover, the operator determined by (|3.15p corresponds 
the operator determined by relation (|3.7D 

A nl (x,t) = A nl (x + X u (t),t), 
which follows from M 3 (t, s)(X 1A - X 7 ) = -AM 3 (t, s)(X 1A - X 7 ). 

4 Conclusion 

Symmetry analysis of an equation is usually performed when solutions of the equation are not 
known, and the basic purpose consists in finding as wide as possible classes of partial solutions 
by using the symmetries of the equation. 

It should be noted that a direct calculation of symmetry operators for a nonlinear equation is, 
as a rule, difficult because of the complexity of the determining equations [32]. The basic subject 
for study is the symmetries related to the generators of one-parameter subgroups of a Lie group 
of symmetry operators [3]. The determining equations for the symmetries are linear, and to solve 
them, it is necessary to set the structure of the symmetries. Finding of nonlocal symmetries 
for differential equations or differential symmetries for nonlocal equations from the determining 
equations, faces mathematical problems. 

In this context, the algorithm proposed in this work to calculate the symmetry operators for 
the FPKE in explicit form is of interest as it provides a possibility to consider the properties 
of symmetry operators for a nontrivial nonlinear equation. The algorithm is stated in terms of 
a direct nonlinear evolution operator and the corresponding left-inverse operator, which enables 
one to vary the structure of the obtained symmetry operators. 

As for the considered equation the evolution operator is found, the symmetry operators are 
not of interest for finding of partial solutions of the equation. However, the symmetry operators 
can be used to study the properties of the solutions obtained. Furthermore the algebra of 
symmetry operators of the linear FPKE (|2.9p (see, e.g. [33]) can be used to study the symmetry 
operators of the nonlinear equation ()2.2|) . 

The developed approach permits a generalization for FPKE's with a smooth operator symbol 
of arbitrary form. In this case, we deal with a solution of the FPKE approximate in the small 
parameter e [21] . The FPKE considered in this work arises in constructing the leading term of 
semiclassical asymptotics for a FPKE with arbitrary coefficients. 
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